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The chirally improved (CI) fermion action allows us to obtain results for pion masses down to 
320 MeV on (in lattice units) comparatively small lattices with physical extent of 2.4 fm. We use 
differently smeared quarks sources to build sets of several interpolators. The variational method 
then leads to excellent ground state masses for most mesons and baryons. The excited state signals 
weaken in quality towards smaller quark masses. In particular the excited baryons come out too 
high. 
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1. Study with two dynamical CI fermions 

We are presenting results of a hadron mass spectrum calculation with emphasis on possible 
identification of excited states. The gauge configurations have been obtained with dynamical, mass 
degenerate up and down quarks. For the fermions we used the chirally improved (CI) Dirac operator 
Dqi tfQ This is a parameterized fermion action of the form 

16 

D mn = £ r a £ c a p n Ui 5n, m+p , (l.i) 

a=l pe^«„ lep 

where &® n symbolizes paths from site m to n. Inserting the ansatz in the Ginsparg-Wilson (GW) 
equation, truncating the length of the contributions (to, e.g., distance 4), and comparing the coef- 
ficients, leads to a set of algebraic equations, which can be solved (norm minimization). We used 
a truncated action with terms involving coupling to neighboring sites within a 3 4 hypercube plus 
some extra terms, giving rise to several hundred terms. Part of our definition of the Dirac action 
is one step of stout smearing ^ of the gauge configuration. The eigenvalues of Dqi are closer 
to the unit circle (where the eigenvalues of exact GW-operators are located) than those of, e.g., 
the improved Wilson operator for the same lattice size (in lattice units). However, (lattice) chiral 
symmetry is still violated, albeit to a smaller amount. The gauge action used is a tadpole-improved 
Liischer-Weisz action. 

The dynamics was implemented with the hybrid Monte Carlo (HMC) algorithm with Hasen- 
busch mass preconditioning (with two pseudofermions) and a chronological inverter utilizing the 
mixed precision technique Details and parameters of the action and the methods of simulation 
can be found in [Q]. There one also finds a discussion of equilibration and determination of lattice 
spacing and AWI (PCAC) mass of the quarks. 

The analysis presented here is based on three ensembles of gauge configurations for lattice 
size 16 3 x 32, with parameters summarized in Table [[]. 



set 




aniQ 


tMD 


configs. 


a[fm] 


m n {MeY] 


BiAwtMeV] 


A 


4.70 


-0.050 


600 


100 


0.151(2) 


525(7) 


42.8(4) 


B 


4.65 


-0.060 


1200 


200 


0.150(1) 


470(4) 


34.1(2) 


C 


4.58 


-0.077 


1200 


200 


0.144(1) 


322(5) 


15.3(4) 



Table 1: Overview of the three ensembles of gauge configurations on which this analysis is based. The 
parameters given are: gauge coupling, bare mass parameter a mo, number of molecular dynamics (MD) time 
units, number of analyzed configurations (in equilibrium and each separated by 5 MD units), lattice spacing 
a determined from the static potential with Sommer parameter 0.48 fm, pion mass, AWI mass (from the 
PCAC relation). For further details see 

2. Variational analysis and hadron interpolators 

In the variational method ^ one studies the cross-correlation matrix of several lattice oper- 
ators Oi with the correct quantum numbers. Inserting a complete set of states (and assuming a 
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discrete spectrum, which is always the case for finite lattices) one finds a superposition of expo- 
nentially decaying contributions, 

C(t)ij = (Oi(t)0](p)) = £(0\Oi\n}(n\o)\0}e- tM " . (2.1) 

n 

Assuming that we have a complete enough set of operators the solution of the generalized eigen- 
value problem 

C{t)vi = X i {t)C{t )v i (2.2) 
allows one to disentangle the individual states and the corresponding energies, 

A,-(0 oc e-' M ' (l + ^(e-' AM ')) . (2.3) 

The eigenvectors are "fingerprints" of the states which one may follow through several time slices 
in order to ensure that the state has been identified consistently (see also the discussion in [[7|].) 

In order to increase the number of hadron operators and to improve the correlation signal 
quality we built interpolating fields with different smearing of the quark fields. We used Jacobi 
smeared quark sources, e.g., u s = S s u, with an hermitian smearing operator S s as discussed in [@], 
with parameters adjusted to produce two different smearing widths, a wide source (s = w, radius 
0.55 fm) and a narrow source (s = n, radius 0.27 fm). 

We also used derivative quark sources as discussed in e.g., sources like ug = DkS w u, 
where k denotes the spatial direction of the covariant derivative 

Di (x, y) = U t (x, 0)S(x + i, y) - U t (x - i, 0) 1 8 (x - i, y) . (2.4) 

The interpolating field operators are built on (in the 3D time slices) 3 times HYP smeared 
gauge configurations []|] with smeared valence quark sources u s , d s (and the strange quark s s ). We 
regularly shift the center of the sources when passing from one configuration to the next in order to 
improve decorrelation. 

For the meson operators we use bilinears. Depending on the quantum numbers this allows for 
sets of operators with different Dirac structure and varying spatial extent. As an example, for the 
pseudoscalar meson we have 

u„Y 5 d n , u n Y5<l w , TiwY5d w , u n y t j 5 d n , u n YtY5d w , u w y t y^d w , 
TiYs d n , u d . YY5 d w , u d , YiYYs d n , u d . YiYtYs d w , ug, 75 d d . , u d . YtYs dd, ■ (2-5) 

The nucleon interpolators have the form 

n® = £abc rf u a (u T b rf d c - 4 rf Mc ) (2.6) 

with the choices (T®, ) = (1, C75), (75, C), and (il, C7475) for i = 1,2,3 respectively, where u, 
d denote again smeared quarks. The A interpolator is 

O^k = £abc u a (u T h CYkU c ) , k = 1, 2, 3 . (2.7) 
projected to spin |. The baryon propagators are projected to definite parity. 
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Figure 1: L.h.s.: The ground state pion and its first excitation. Open symbols denote partially quenched 
data, the full symbols denote the dynamical data for the three ensembles A, B and C. R.h.s.: The dynamical 
data for the ++ channel ao are compared with the energy level expected for a % — r\2 channel. In the insert 
also the partially quenched data (open symbols) are shown, as discussed in the main text. 



3. Results for mesons 

+ : 7T ± (140), 71^(1300). As discussed in [BJ|, in the multi-operator (variational) analysis at 
small pion masses the backwards running (in time) pion limits the observation range for the excited 
state. This can be cured by a larger time-size; however, for physical pion masses we expect that 
one needs at least N t = 64 for lattice spacing a = 0.15 fm or N t = 128 for a = 0.075 fm. 

We choose the fit interval of the exponential fit to the non-leading eigenvalues based on the 
window where the backwards running contribution is not yet dominant and get the excited pion 
signal in Fig. [j] (left). We include in this plot (as in some of the other figures) the results for 
partially quenched data, i.e., where the valence quark masses are larger than the sea quark masses. 

++ : ao(980), ao(1450). The isovector, scalar meson has led to controversial results in lattice 
simulations JTo| , [TTJ] . Most quenched studies found a ground state extrapolating towards the mass 
of the £?o(1450) for smaller valence quark masses. Results for (two) dynamical quarks seem to 
lead to smaller masses, compatible with an extrapolation towards the £?o(980). However, for these 
masses the energy values are close to those of an expected % — f]2 channel in s-wave (mass of 772 
estimated [|I2|]). Fig. [T] (r.h.s.) exhibits the situation. 

When also plotting partially quenched values, we find an interesting effect for the ensemble C 
with the smallest sea quark mass. The partially quenched data do not smoothly extrapolate to the 



dynamical point. An explanation has been offered in [13]: the partially quenched states may couple 
to pairs of pseudoscalars (composed of valence and sea quarks), leading to unphysical contributions 
that cancel in the fully dynamical case. 

Also, we find a broad range of values for the extracted energy levels for different sets of 
contributing interpolators. The issue, whether the ao ground state is dominantly a tetraquark state, 



is still not settled (see [14]) 



1 : p ± (770), p ± (1450). In Fig. |2| (l.h.s.) we show the results for ground state and first ex- 
citation in the p -channel. For a better signal we allow for different operator combinations for the 
excited states. The broken lines give an error window for linear extrapolating fits based on the three 
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Figure 2: L.h.s.: Ground and excited state in the 1 channel. Open symbols denote partially quenched 
results. R.h.s.: Ground state in the 2 ++ channel separately evaluated for interpolator sets in the E (top) 
and the representation (below) of the cubic group. Filled symbols denote the values from the three fully 
dynamical data sets, open symbols show partially quenched results. 

dynamical points only. The agreement with the experimental values is surprisingly good in view 
of the p being a resonance; however, the decay pions are in relative p-wave and due to the given 
lattice size the necessary extra unit of momentum stabilizes the vector meson. 

We find more admixtures of still higher excitations, and for some combinations of operators 
we can identify a 3rd energy level compatible with the p(1700). 

Our data in the exotic 1~ + channel is too noisy to allow extraction of a ground state energy 
level. 

1 ++ : ai(1260) and l + ~ : bi(1235). In both channels we see reliable signals only when includ- 
ing derivative interpolators (cf. the discussion for the quenched case in [^]). The error bars are 
somewhat larger for l + ~ than for 1 ++ , where they are roughly ±50 MeV. In both cases a lin- 
ear extrapolation points towards the experimental value (within the errors). The mass difference 
miat,b { ~ miatt,a Q is approximately 200 MeV. 

2 ++ : a2(1320). In this channel two representations of the cubic group couple: E and T%. Again 
interpolators involving derivative sources are necessary, e.g., EjjkUg.jjdn. We find good signals 
only when we include several interpolators and thus this is a case where the variational method is 
crucial. Fig. ^| (r.h.s.) demonstrates the situation. The results for the independent analysis of both 
representations are compatible, although T2 appears to extrapolate better to the physical value. 

2~+ : 7T 2 (1670) and 2 : p 2 (1940)?. Qualitatively the situation is like for the 2++ channel. In 
both representations we find values compatible with each other and within large errors of order 200 
MeV also (in linear extrapolation) with the physical mass value. 

4. Results for baryons 

| + : N(940), N(1440), N(1710). The baryon system poses various challenges in the excited sec- 
tor, among them the identification of the Roper state in lattice calculations. In most calculations in 
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Figure 3: Results for the positive parity nucleon channel for the new data with dynamical fermions (l.h.s.) 
compared to earlier results from quenched configurations JT^ ] (r.h.s.). 
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Figure 4: Positive parity: Ground state A and first excitation. 



the quenched approximation the first excitation of the (positive parity) nucleon came out too high, 
and only within large error bars an extrapolation towards the experimental Roper mass was imag- 
inable (see, e.g., the quenched study in [|I5|]). Only recently, and at this conference Jl6| , 17], there 



are results showing an energy level extrapolating towards smaller masses. That study was based on 



optimizing combinations of interpolators with several widths, like in this and earlier work Q18p. 

There was some hope that introducing the quark dynamics might improve the situation. Our 
present analysis does not strengthen that hope (Fig. ||). Indeed our new results for dynamical 
fermions are similar to earlier quenched results [[HI]. We want to emphasize that it is important 
in this case to simultaneously find both excitations (the Roper and the state extrapolating to the 
A^(1710)) in order to have a convincing identification. Both excited states are too high in our 
results which might indicate that the given volume is too small for baryon excitations. 

2 : N(1535), N(1650). We clearly identify two states but cannot quantify their mass splitting 
due to large errors. 

| + :A(1232),A(1600) and §~ :A(1700). In the positive parity sector (see Fig. g) we find that 
the ground state is closer to the experimental values, as compared to the quenched results. The 
excited state is clearly seen, but again too high, maybe due to volume squeezing. The ground state 
in the negative parity sector is clearly identified and extrapolates towards the experimental mass. 
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Thus, while our meson results look quite good, the excited baryons come out too high. Im- 
provement of the signal might be obtained with larger lattices and possibly extending the set of 



interpolator (see, e.g., [|19|]). Further results for the axial charge of baryons including baryons with 



strangeness are discussed in []2Q|]. 
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